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In this paper we present a global fit of beyond the Standard Model (BSM) dimension six operators 
relevant to the top quark sector to all currently available top production cross-section measurements, 
namely parton-level top-pair and single top production at the LHC and the Tevatron. Higher order 
QCD corrections are modelled using differential and global A'-factors, and we use novel fast-fitting 
techniques developed in the context of Monte Carlo event generator tuning to perform the fit. 

This allows us to provide new, fully correlated and model-independent bounds on new physics 
effects in the top sector from the most current direct hadron-collider measurements in light of the 
involved theoretical and experimental systematics. As a by-product, our analysis constitutes a 
proof-of-principle that fast fitting of theory to data is possible in the top quark sector, and paves 
the way for a more detailed analysis including top quark decays, detector corrections and precision 
observables. 


I. INTRODUCTION 


The Standard Model of particle physics has proven to 
be an extremely successful description of Nature up to the 
electroweak scale. Nonetheless there are many compelling 
reasons to believe it is an intermediate step to a more 
fundamental picture of physics at the TeV scale. 

The top quark, as the heaviest Standard Model particle, 
is expected to play a unique role in this new physics. Given 
the unsatisfactory explanation of electroweak symmetry 
breaking within the SM and the appearance of rrit at the 
electroweak scale, i.e. the closeness of the top Yukawa 
coupling to unity, the top mass may arguably be seen as 
a strong hint of physics beyond the SM. 

Most BSM scenarios lend a special role to the top quark. 
In supersymmetry the light Higgs mass is stabilised from 
UV divergences by the contribution of SUSY top partners, 
among others (see, e.g.. Refs. mm)- In compositeness 
scenarios mm, the quark masses and EWSB are gener¬ 
ated through linear couplings of the SM fermions to new 
strongly-interacting physics at the TeV scale. In theo¬ 
ries of warped extra-dimensions, the top quark couples 
preferentially to Kaluza-Klein states in the 5D bulk [5JI1], 
offering a unique window to the new physics. 

Typically all these scenarios predict a modification of 
Higgs phenomenology, which has been thoroughly stud¬ 
ied after the Higgs discovery IZHII]- Such analyses are 
currently limited by small statistics in the observed Higgs 
discovery modes. Taking the special role of the top quark 
in electroweak symmetry breaking at face value, the abun¬ 
dant production of top quarks at the LHC provides a 
complementary avenue to search for new non-resonant 
physics beyond the SM, that will be relevant to our un¬ 
derstanding of electroweak symmetry breaking. 

Given the plethora of concrete scenarios and the ab¬ 
sence of any telling signals of new physics in the current 
data, parametrising BSM effects in an effective field the¬ 
ory expansion [28j is well-motivated. In this approach, 
all possible interactions are captured in an effective La- 


grangian Ces: 

•^eff = -^SM + ■■■ 

The higher-dimensional Lagrangian terms Ci are sup¬ 
pressed by powers of A - the energy scale associated with 
the new physics. In the top-down approach, we have 
integrated out all heavy degrees of freedom, capturing 
their low energy phenomenology guided by SM gauge 
and global symmetries, irrespective of their concrete UV 
dynamics. Such an expansion is valid provided there is a 
good separation of scales between the typical collider en¬ 
ergy and A. However, this approach is completely general: 
the {Ci\ are constructed from SM operators, respecting 
the 51/(3) X SU{2) x [/(I) gauge symmetry. 

The leading contributions relevant to new physics in 
the top sector enter at the dimension-six level Oil/h?) 

Aff = -CSM + ^ X ’ 

i 

where Ci are arbitrary ‘Wilson coefficients’ and Oi are 
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TABLE 1: All dimension-six operators relevant to top quark 
production, in the notation of Ref. m- Details of each are 
included in the text, q denotes the left-handed quark doublet, u 
and d denote the up-type and down-type right-handed singlets. 
We do not include explicit flavor indices here, the relevant 
flavor indices are included in the text. 13 operators are shown, 
but Otw and OtG have both real and imaginary parts which 
should be considered as independent operators; the latter 
produce GP-violating effects. 
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Dataset 

y/s (TeV) Measurements 

Ref. 

Top pair production 



ATLAS 

7 -b 8 

Total inclusive a 

m 


7 -b 8 

Differential Prit), Mu, \y{tt)\ 

m 

CMS 

7 

Differential pT(t), Mu,y(t), \y{tt)\ [1^ 

CDF 

1.96 

Differential Mu 

m 

D0 

1.96 

Differential Mu,Prit), \y{t)\ 

m 

Single top production 
ATLAS t-channel 7 

Total inclusive a 

[m 


7 

Differential prit), \y{t)\ 

CMS t-channel 

7 

Total inclusive a 

m 


8 

Total inclusive a 

m 

CDF s-channel 

1.96 

Total inclusive a 

m 

D0 s + t-channel 

1.96 

Total inclusive a 

m 


TABLE II: Datasets used in the fit, including 
total cross-sections (a); transverse momenta 
of single tops (pt(I)) and top pairs {pT{tt))\ 
rapidities of single tops {y{t)) and top pairs 
(y{ti))-, and the invariant mass of top pairs 
{Mu). 


dimension-six operators. These operators lead to no¬ 
ticeable deviations from SM expectations in a double 
expansion of the matrix element in SM and new physics 
couplings 

lA^totP = |A^SMp + 23?{AIsMA^r)6} + (1) 

where strictly speaking one must neglect the third term 
on the right-hand side if working to dimension-six only, as 
this has dimension-eight, . Provided C^/A^ is small, such 
a truncation is typically valid and the squared dimension- 
six terms become numerically irrelevant. 

The complete set of 80 effective operators at dimension- 
six has been known for some time Only recently 

was it shown that this basis contains several redundancies, 
with the minimal set comprising 59 terms [HHZlllHlllH]- 
Considerable attention has been devoted to constraining 
these operators, for example, in the context of Higgs and 
precision electroweak physics [IHII]. In addition, strong 
bounds have also been placed on new top interactions 
from precision constraints at LEP [5S] and direct searches 
for top quark physics at the LHC [30U35] . 

While Higgs physics has received a lot of attention 
from an EFT perspective, the top quark sector has not 
seen similar scrutiny, although top data from the combi¬ 
nation of the Tevatron and the LHC Run I is far more 
abundant. In the last few years, top quark physics has 
entered something of a precision era: the top has been 
measured in several production and decay channels, and 
dedicated searches in complicated final states such as tiH 
are underway [ 311371 . 

It is our aim to close this gap. The TopFitter ap¬ 
proach constrains new physics in the top sector using 
both differential and inclusive observables, by means of 
a computational tool which is fully flexible with respect 
to the number of input measurements and scales well to 
the relevant number of EFT operators. In the present 
work we limit ourselves to a nine-dimensional fit based on 
direct top measurements performed at the Tevatron and 
the LHC, keeping track of all EFT operator-correlations. 


and reserve a more complete investigation for the near 
future [38] . 


II. RELEVANT OPERATORS 

Throughout the analysis, and for ease of compari¬ 
son with precision electroweak studies, the operator set 
presented in Ref. [12] is used (see also the basis of 
Refs. [331l3ni)- Assuming minimal flavor violation, and in 
the leading-order* approximation of equation 2, of these 
59 operators only 15 — shown in Table [I]— are relevant 
for top production. Fitting a 15-dimensional function is 
a considerable challenge; a brute force likelihood scan at 
N points per dimension would require evaluations, 
which is prohibitive even for modest, low-resolution values 
of N. This naive dimensionality can be reduced, however, 
by noting some features of the operator set. 

Firstly, we note that the two operators containing the 
dual field-strength tensor , along with 

the imaginary parts of Otc and OtWi are CV-odd and 
can be discriminated from CT^-even effects in studies of 
spin-correlations, polarisation effects and genuinely CV- 
sensitive observables [42] (for recent analyses focusing on 
the tWb vertex, for instance, see Refs. [43H47] i. Currently 
there is no evidence for CT^-violation in the top sector 
beyond the minimal flavor violation assumption. We will 
address these operators in forthcoming work but neglect 
them in the following; the dimensionality of our fit is 
reduced by four. 

Secondly, we consider top-pair production. Here the 
four-fermion operators, which are numerous when all 
flavour combinations are considered, only contribute to 


*By leading-order we mean (P(A“^), but for some new physics 
effects, such as top flavour-changing neutral currents, the first non¬ 
zero contributions enter at 0(A“'^) see e.g. m for details. 
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top-pair production through the partonic subprocesses 
uu, dd —>■ tt, which reduces the myriad of possible op¬ 
erators to four unique, flavour-specific linear combina¬ 
tions 15^1 H51: 


ill33 

'‘uu 


c; = 3(2c(«‘“ + c;”) 

_ (Ml) 1133 I ^(3) 1133 I f-iL 
\^qq ^qq ' 

C3 = 3(cg>™‘-cm‘”) 

+ (cg>‘'»»-c<;>‘'») + 6c™”“ 

^2 _ /^(8) 1133 , ,-.(8)3311n 

'-'d — ~ I'-'nti + '■^nd b 


where explicit flavor indices q^) have now been 

included. The non-4-fermion operators Og, Otc^ and 
O^G also contribute to top pair production, giving a 
total of 7 relevant operators. In the gg —>■ tt channel, 
Og rescales the triple gluon vertex while OtG modifies 
the top-gluon coupling; only contributes through 
gg ^ h ^ tt, which is heavily suppressed in the Standard 
Model although it can be probed in ttH production. 

Three CP-even operators^ contribute to single top pro¬ 
duction: Otw modifies the tWb vertex, as does while 

the operator creates a new four-quark topology 

which interferes with the SM piece. 

There is hence a clean factorisation into 7-1-3 CP-even 
operators associated with top quark production at hadron 
colliders. In this study we reduce this further to a 6-1-3 
configuration by eliding the highly suppressed contribu¬ 
tions of 0(f,G to top-pair production. 


III. DATASETS 

The aim of this paper is to present a preliminary study 
demonstrating the feasibility of performing a full global 
fit of top quark effective theory to data. We thus include 
top quark pair and single top production processes at 
parton level only, which observables and datasets 
are collected in Table UTl 

It should be noted that a fully differential fit along these 
lines consists of a multitude of exclusive measurements. 
Treating each bin as an independent^ measurement, we 
have 103 bins for top pair production and 23 from sin¬ 
gle top. This highlights the necessity of a fast analysis 
framework, as introduced in the present work. 

Given that we will model higher-order corrections as 
described in the following section, we do not include Wt 
production, which interferes with top pair production 


tThe contribution of the operator is 

heavily suppressed, as its interference with the SM amplitude is 
proportional to mi, (see e.g. Hi)- 

tWhere published by the experiments, we have included bin-to- 
bin correlations. These have negligible effect on our conclusions. 


at next-to leading order, such that it is not possible to 
reproduce existing experimental analyses using a fixed 
order parton level calculation [SOHMj . 


IV. DETAILS OF ANALYSIS 

We begin by including the operators listed above (to¬ 
gether with consequent SM parameter redefinitions) in a 
FeynRules m model file, which is then interfaced via 
UFO [56] to MadGraph/MadEvent |56l[57] in order 
to obtain parton-level theory predictions. Samples were 
generated for all the relevant processes: top pair produc¬ 
tion: pp tt, single top production: pp tb (s-channel), 
and pp —>■ tq (t-channel). 

In order to model next-to leading order QCD correc¬ 
tions, SM-only samples at next-to-leading order are gen¬ 
erated with Mcfm [^. These are used to construct 
differential (bin-by-bin) and global A-factors, as in e.g. 
Ref. [59]. Theoretical uncertainties for these samples are 
estimated in the usual way, by independently varying 

the scales I^central/^ ^ pK.F ^ ^Mcentrab where /^central 

is taken to be m*. Parton distribution function (PDF) 
uncertainties are estimated by generating events using the 
next-to-leading order NNPDF23 (5D], MSTW2008 [5T] . 
and CTIO |S5] PDF sets, according to the PDF4LHC [55] 
prescription. We take the central value as our estimate 
and the width of the envelope (including scale variations) 
as the total theoretical uncertainty. In the case of top 
pair total inclusive cross-sections, we use global iF-factors 
from next-to-next-to leading order QCD with soft gluons 
resummed to next-to-next-to-leading logarithmic accu¬ 
racy [5iH55] . 

A strength of our fitting procedure is the use of novel 
techniques developed in the context of Monte Carlo event 
generator tuning, as implemented in the PROFESSOR [55] 
framework. The procedure is as follows: 

• A set of points in the V-dimensional parameter 
space {Ci} is sampled logarithmically. Other sam¬ 
plings are possible - we choose logarithmic sampling 
to avoid oversampling of regions where coefficients 
are large, such that dimension-eight terms become 
important. 

• At each sampled parameter space point, all the¬ 
ory observables are calculated, with uncertainties, 
as described above. One then constructs a poly¬ 
nomial parameterising function fbi{Ci}) for each 
observable bin b, which fits the sampled points with 
least-squares-optimal precision. This function can 
be used to efficiently generate theory predictions for 
arbitrary parameter space points within the fitted 
range. We choose a third-degree polynomial for 
this function. This has been shown to work well 
in Monte Garlo tuning [69], and should in fact be 
better suited to the present case: in the absence 
of uncertainties, each observable is a second-order 
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individual 

Cg i-^—i marginalized 

CtG , _ 

Cl \ , __, 

cl \ ' 

, 

m ^ ^' : • '^ 

Ctrv , ' , 

C^q “h 

Clq ~ , 


-1 -0.5 0 0.5 1 

Ci = 

FIG. 1: 95% confidence intervals for operators contributing 
to top-pair and single top production, individually (with all 
other operators set to zero) and marginalised (with all other 
operators allowed to float to best-fit values). Note that the 
marginalised bound on Cj fall outside the region where the 
dimension-six approximation is valid, so this operator is un¬ 
constrained. 


polynomial in the {Ci}, cf. eq. Q, and the ex¬ 
tra polynomial order provides some tolerance to 
beyond-fixed-order effects. 

• Finally, we construct a function between the bin 
parameterisations {fbi{Ci})} and the data, accord¬ 
ing to 


O b 


(A({CJ) - 
^2 


i.e. we sum over all observables O, and all bins in 
that observable, b. Eb is the experimental reference 
value at bin b and ab is the total uncertainty for bin 
b, which we for now assume as an uncorrelated com¬ 
bination of theoretical modelling and experimental 

measurement uncertainties, ab = o'^ieory + ^exp■ 

The is then used to place constraints on the 
operator Wilson coefficients, as follows. 


Constraints are obtained in two ways, for ease of com¬ 
parison with existing literature. Firstly, single operator 
coefficients are allowed to vary, with all others set to 
zero (the SM value). The x^ is then minimised using 
PyMinuit m , and used to set confidence limits on the 
operator value. A second approach is to marginalise over 
the remaining operators, namely to construct the confi¬ 
dence limit for a given operator coefficient whilst allowing 
all other coefficients to vary. Both cases are shown in 
Figure where the dimension-six contributions are nor¬ 
malised to the Standard Model piece via Ci = CiV^ jIs?. 
All results are consistent with the SM within 95% limits. 


As with all effective operator constraints, these must 
be interpreted as valid only in the region where 0{hr'^) 
terms are not large. Clearly Cj is outside this region. In 
top pair production, for instance, the contribution from 
dimension-6 operators relative to the SM piece is typically 
0{glCiv ‘^!which must be < 1 in the linear approx¬ 
imation , i.e Ci < 1.5. All other operators respect this 
bound. It should be noted that some of these operators, 
namely those containing field strength tensors, can only 
be generated at loop level in the ultraviolet completion, 
which widens this region of validity since A^ will be ac¬ 
companied by a loop factor of 167r^. This argument is 
invalid, however, if the underlying completion is strongly 
coupled. It is possible to include such information in our 
fitting approach, but in the interests of full generality no 
such model-specific assumptions are made here. 

One sees from Figure that the weakest constraints 
are on the coefficients (of four-fermion operators) C^ and 
(7^. These are constrained by the processes uu —>■ ti and 
dd ti respectively, which are suppressed relative to the 
corresponding gluon initiated processes, mostly due to 
the relative partonic luminosities. 

One may also examine the correlation of constraints 
between pairs of operators. An example is Figure |2(a)[ 
which shows confidence limits in the {CtWiC^^) plane, 
with all other operator coefficients set to zero. One may 
also marginalise over all remaining operators, as shown 
in Figure [2(b)[ In both cases, we currently find excellent 
agreement with the SM. More detailed results will be 
presented in a forthcoming paper [35]. 


V. SUMMARY, CONCLUSIONS AND 
OUTLOOK 

Following the discovery of the Higgs boson, the search 
for physics beyond the Standard Model will remain the 
primary goal of the LHC experiment for the foreseeable fu¬ 
ture. The top quark sector is a particularly well-motivated 
window through which to look for the imprint of non¬ 
resonant new physics. Modelling such effects using effec¬ 
tive field theory (higher dimensional operators) is well 
justified given the absence of new resonant physics from 
the LHC Run 1. The abundance of top quark production 
at the LHC enables a multi-faceted analysis of top quark 
phenomenology and allows us to confront higher dimen¬ 
sional top sector operators with differential measurements 
at high statistics. 

In this paper, we have characterised new physics cor¬ 
rections using the well-established framework of effective 
field theory. We have presented results from a new com¬ 
putational framework to fit all possible dimension-six 
operator coefficients to a comprehensive set of relevant 
data. This is possible through our use of fast-fitting algo¬ 
rithms, which have been developed (and well-tested) in 
the context of Monte Carlo event generator tuning. Here 
we expect these techniques to work even better, given the 
explicit polynomial dependence of theory observables on 
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(a) individual constraints. (b) marginalised constraints. 

FIG. 2: 68% (blue), 95% (turquoise) and 99% (orange) confidence intervals for Ctw and C^q in a global fit, with all remaining 
coefficients set to zero (a) and marginalised over (b). The star marks the best fit point, indicating a currently good agreement 
with the Standard Model. 


operator coefficients. 

Our method involves constructing a parameterising 
function to effectively parametrise the theory output of 
Monte Carlo generators (here at parton level only). Once 
this has been constructed, it is quick to perform a global 
fit containing all possible operators, and to amend this fit 
as and when new data appear. Furthermore, there is no 
significant speed decrease in our fitting procedure upon 
improving the theory prediction (e.g. to include parton 
shower or detector corrections), as such improvements 
only affect the parameterising function, which has to be 
calculated only once. 

The results of our fit currently show good agreement 
with the Standard Model, which is unsurprising given 
the absence of new physics currently reported in other 
studies. Our results, however, provide a proof of princi¬ 


ple study that efficient global fits of top quark effective 
theory are possible. It is straightforward to generalise 
our fit to include more experimental observables (beyond 
parton level, including top quark decays), to improve the 
theory description with higher order corrections, and to 
include new data sets including those from the recently 
commenced LHC Run II. Work in these directions is 
ongoing. 
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